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Abstract. - 

We study the anomalous scaling of the mass density measure of Lagrangian tracers in a 
compressible flow realized on the free surface on top of a three dimensional flow. The full two 
dimensional probability distribution of local stretching rates is measured. The intermittency 
exponents which quantify the fluctuations of the mass measure of tracers at small scales are 
calculated from the large deviation form of stretching rate fluctuations. The results indicate 
the existence of a critical exponent Uc — 0.86 above which exponents saturate, in agreement 
with what has been predicted by an analytically solvable model. Direct evaluation of the multi- 
fractal dimensions by reconstructing the coarse-grained particle density supports the results for 
low order moments. 



Advection of homogeneous distribution of passive particles in a compressible flow gener- 
ically results in clusters of particles' concentration [1-3]. The effects of compressibility in 
three-dimensional flows are typically relevant only at Mach numbers near or larger than 1 
when density inhomogeneity reacts back on the fluid velocity. However there are simple 
physical systems in which the phenomenology of passive particle advection in compressible 
flows is applicable. Effective flows of incrtial particles in incompressible flows at small Stokes 
numbers [4-6] and flows of surface suspensions [7-9] are among such instances. 

Consider the distribution of particles advected by a velocity field u(x, t) assumed homoge- 
neous and stationary in time. The evolution of particle density p(x, t) in the limit of vanishing 
diffusivity is given by 

§ + V-(pu) = (1) 

If V • u ^ the flow is compressible and particle trajectories asymptotically converge to a 
non-uniform distribution. Compressibility can be characterized by the dimensionless ratio 
C = iif^ which assumes values in < C < 1, the limiting values corresponding to an in- 
compressible and a potential flow respectively. The particle density depends on the particular 
realization of the velocity field and for a stationary random flow the density is expected to 
converge to a statistically stationary state p*. 
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A suitable way to characterize the statistics of the density field on the attractor is to 
measure the Lagrangian moments of the coarse-grained mass in a small region of size r as: 



(2) 



where mr(x) — p*(x')(ix' is the coarse-grained mass in a ball of radius r centered at 

X. The weight in Q restricts the average to balls centered on the attractor. Note that 
the coarse grained mass mr(x) depends on the specific velocity configuration within a given 
ensemble of flow realizations: the (• • ■) means the averaging with respect to the given velocity 
statistics. 

When r — > and under generic assumptions for the random flow (m") ~ r^" [3, 10]. A 
non-linear dependence of the set of scaling exponents ^„ on the order n implies an intermittent 
distribution of mass fluctuations at small scales. Recently, a phenomenological connection be- 
tween the probability distribution of small scale velocity stretching rates and scaling exponents 
of coarse grained mass moments has been obtained [11]. More precisely the set of exponents 
^„ are bridged to the Cramer function of the stretching rates of the underlying flow. It is 
shown that ^„ saturates at large n. The saturation is a footprint of the very violent events 
which spoil the scaling in the higher mass moments. 

Here we study the clustering of passive tracers which move on the two-dimensional free slip 
surface on the top of a three-dimensional incompressible fluid [8,9]. Our main result is that an 
intermittent distribution of passive tracers, with saturation of high-order scaling exponents, 
is observed also in this realistic, time correlated flow. 

When the flow is spatially smooth the evolution of trajectories is a combination of random 
stretching and contracting in time. The individual tracers are advected by the same realization 
of the flow thus the evolution of the separation between two Lagrangian particles is given 
by X(t, x) ~ W(i,x)X(0,x) where X(0,x) is the initial separation. The linear random 
stretching and contracting is encoded in the positive eigenvalues of the symmetric matrix 

{t, x)W{t, x): e^*'^i , • • • , The forward stretching rates ai{t) are arranged in non- 

decreasing order and asymptotically converge to the Lyapunov exponents at large times, i.e. 
Ai = limt^oo (7i [3]. For generic cases when the spectrum is non-degenerate the local stretching 
rates satisfy an asymptotic probability distribution function (PDF) with the following large 
deviation shape [12], 

P(ai,...,arf;t)^e^*^(-^'--'^''), (3) 

where the Cramer function H(ai, • • • , ad) > is convex and takes its minimum in H{\i) = 0. 

For completeness, here we briefly sketch the outline of the theoretical arguments first 
appeared in [11]. In a two-dimensional compressible flow, the Lagrangian evolution correspond 
to a dissipative dynamics with (Ti > and cr2 < and with Ui < 0. Along the unstable 
direction the density is smooth while in the stable direction the density conflguration is fractal. 
After a transient regime the density is expected to reach to statistical stationary state. We 
will take formally the initial time to —oo thus we will assume that aX t = the density 
is already stationary. At time t > one looks at the fluctuation of the coarse-grained mass 
generated by p, according to Consider a rectangular box of size r with sides ri and r2 

£(1) £(2) 

along the orthogonal unstable and stable directions. According to eq.® ^2) ~ '^1" ^2 
where ^„ = ^l-P + ^i?''. Since the measure along the unstable direction ri is regular one gets 
= min(^„,n). For ^„ > n it results to ^i^^ — n and therefore (m"^ ^ r"r2"~". In the 
case of ^„ < n it reads {m"^ ,.^) ~ rfv§. 
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The pre- image of the box at time i = is a stretched box as ri shrinks to rie^*"^ and r2 
increases to r^*"^^ . Due to conservation of the mass and stationarity the mass moments in the 
pre-image box at time t = are the same, i.e. 



'1 ,r2e 



(4) 



According to the definition in ^ the average on the right hand side is considered over the 
distribution of the forward stretching rate with respect to invariant measure p* . Note that the 
result of averaging over forward and backward distributions with respect to invariant measure 
are the same [13]. At large time t one resorts on the saddle point argument whereby two 
distinct situations emerges from eq.Q that is, 



mm 



in {i/(cri, 0-2) + "CTi + - 71)0-2} = n<^n, (5) 



min{i/((Ti,o-2) + n > (6) 



CTl ,0"2 



The above relations connect the scaling exponents ^„ to the Cramer function H{ai, (72 )■ For a 
synthetic Gaussian and white in time smooth flows (the Kraichnan model of turbulence [14]) 
where the exact functional form of Cramer function is known, eqs. (|5I6|I lead to an analytic 
form for the scaling exponents [11]: 



2n+V(l+2C)^-8Cn^ _ if „ < „ 

_ I 1+2C ^ 11 /t lie , ^^-^ 

Coo if n > Tie 



where the critical order for saturation is given by 



l^J^ + lc if o<c<i, 
if i < C < i 

1+2C 6 — ^ ^ 2 ■ 

In the general case, depending on the flow compressibility C two scenarios for saturation 
are thus conceivable. In the first scenario the line C„ = ^oo from © intersects with the curve 
of Cn('T-) given by ©. In such a case ^^o = Cn^ = and the value of {aKric), 172 ("-c)) will pass 
through the line ai = 172 . In the second case the line tangent to the maximum of the manifold 
C„(n) given by Q identifies ^oo- In this later scenario ^ ric and ric identifies with the 
solution of ^Iric — 0- For the Kraichnan model the second case results in Coo — Cric — 2nc — 1- 

We have integrated the three-dimensional Navier-Stokes equation for an incompressible 
flow by means of a standard pseudo-spectral code at resolution 64^. Lateral boundary con- 
ditions for X and y directions are periodic, while in the vertical z direction we apply free-slip 
boundary conditions, Uz = and dzU^ = dzUy = [7]. When the flow is in statistically 
stationary conditions, 10^ passive tracers are placed on the top free surface z — 0. Their two 
dimensional positions x(t) = (x{t),y{t)) are advected by the two-dimensional flow {ux,Uy) 
which is compressible since dxUx + dyUy — —dzUz ^ [8,9]. The compressibility of the surface 
flow was found to be C ~ 0.21. 

In order to evaluate the stretching rates, we have integrated along each trajectory a set of 
two orthogonal tangent vectors z'^'"') {k — 1, 2) evolving according to 

±zl'^^iV,u.)zf (z,j = l,2) (9) 
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Fig. 1 - Marginal Cramer functions Hi{ai) (left) and 1^2 (0-2) (right) computed at three different 
times lags r — , t = 2A^^ and r — 4A^^. 



The PDF of stretching rates at fixed time lag r, P{ai , cr2 ; r) , is then accumulated by computing 
CTfc = (1/r) ln(|z('^)(T)|/|z('^)(0)|) for different reahzations along the trajectories. 

Stretching rate statistics are accumulated up to t^ax = 15000Aj~^ which correspond to 
about 7000 velocity correlation times T ~ 2.1A^^. The fluctuations of each of the stretch- 
ing rates is quantified by the convergence of their corresponding marginal Cramer functions 
Hi{(Ji). These are defined as 

F,(a,;i) = j P{au'J2;t)da3-^ ~ e-*^'("') ; i = 1,2, (10) 

FigH shows the two marginal Cramer functions obtained from the marginal probability 
density functions of stretching rates at three different time lags t = X'^^ , t — 2A^^ and 
T = 4A]"^ The average values of stretching rates, i.e. the minimum of the marginal Cramer 
functions, converge to the Lyapunov exponents Ai ~ 0.29, A2 — —0.86 from which one es- 
timates the Lyapunov dimension [3] as Dl = 1-1- — 1.34. This value is larger to what 
obtained in a similar set of simulations at higher Reynolds numbers [9] , but is consistent with 
the fact that here the value of compressibility is smaller. 

In Fig.|(2Jl we present the contour plot of the two dimensional Cramer function H(ai,(72) 
evaluated at three different time lags. The overlap of the contour corresponding to different 
time delays is the indication of the convergence of the statistics. We observe that the contours 
are rather far from an elliptic shape, a direct indication of the deviations of the probability 
density P{ai,a2',T) from a Gaussian distribution (as it is also evident from the marginal 
distributions in Fig. P). 

From the Cramer function H{ai,(72) the set of scaling exponents ^„ is obtained by nu- 
merical minimization of relations (jSJ-lO- For each order n this procedure determines a point 
{(7l{n),a2{n)) in the stretching rate plane. The line in Fig.© represents the trajectory of 
these points starting from n = to n = ric — 0.86 where the point crosses the line cti = <72 
and the exponents saturates, for n > ric, to — nc following the second of the two scenarios 
described above. 

The set of scaling exponents obtained from the minimization of Cramer function is shown 
in Fig.|3| Intermittency of mass distribution is evident from the non-linear behavior of 

The geometrical interpretation of the saturation predicted by © is that the intersection 
between the line — ^ricCi s-^id the curve H{ai, a2{nc)) is tangent at n = 71^. In the inset of 
Fig. Owe show that this is indeed reproduced by our data and that the tangent point is well 
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Fig. 2 - Two dimensional contour plots of the Cramer functions H{ai, CJ2) computed at three different 
times r = A^^, r = 2A^^ and r = 4A^^. The geometrical locus of {al{n), (T2(n)) is shown from n = 
(A) to n = ric ~ 0.86 (B). 

inside the convergence range. This resuh indicates that the saturation observed in Fig. is a 
genuine phenomenon and it is not induced by finite size statistics. 

We have compared the scahng exponents determined by the stretching rates with those 
directly measured from the mass distribution. To this aim, we stored 12000 configurations 
of the 10^ particles at an interval At = 0.67T. Particle density p{x) is reconstructed by 
distributing the positions of the Lagrangian tracers on a regular grid of resolution 512 x 
512. For each snapshot we directly measure the coarse-grained mass and calculate the mass 
moments (m") according to Eq.ljSJ. The moments are then averaged over the 12000 snapshots 
and the scaling exponents are estimated by a best fit. The box-counting scaling exponents 
S}n^^^ are shown in Fig.|2|up to order n — 1.3. 

As it is evident from the figure the moments are extremely fluctuating. The low order 
exponents, for n < Uc, are very close to the values obtained from minimization of Cramer 
function. The saturation of the scaling exponents is corresponding to the rare events when 
large fraction of the particles accumulate at very small scales. For n > ric the box counting 
exponents deviate from the saturation line, but we have found that the distribution strongly 
fluctuates from one realization to another. 

In order to estimate the magnitude of these fluctuations, we have computed, for each 
realization, the maximum fraction of particles in a single grid cell, i.e. Hmax/-^- The resulting 
distribution is shown in Fig 0] which shows a rather fat tail with events up to rimax/N ~ 0.4. 
This value of concentration has to be compared with the most probable value Uj^ax/N ~ 0.009 
(see FigEJ. 

These configurations of extreme concentration are responsible for the saturation of scaling 
exponents. Indeed, if we compute the box-counting exponents S,n^'^^ for one of these particular 
configuration, we observe a saturation of the exponents at a value even smaller than (see 
Fig. 131). Although our data provides sufficient statistics to measure the joint distribution of 
local stretching rates the dynamical evolution of the rare events responsible for saturation of 
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Fig. 3 - Scaling exponents computed from the stretching-rate analysis (connected circles). The 
horizontal line for n > Uc — 0.86 represents the saturation value = ric- Scaling exponents obtained 
from a direct box-counting analysis are shown after averaging 60 sets of 200 configurations each (+, 
error bars estimated on the fluctuations among different sets). The lower line is the result of the direct 
computation of mass exponents for the particular subset where saturation is observed (see Fig. 2J ■ 



the coarse grained mass requires by far a longer sampling, 
fluctuation for ^n^^^ > n in Fig. |31 



This explains the large level of 
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Fig. 4 - The distribution of the maximum fraction of particles contained in a single bin of a 256 x 256 
grid. In the inset the time series of riraax/N is shown. The circled point is the "singular" event which 
is contained in the subset from which the dashed curve of Fig. |3 was extracted. 
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In summary we have obtained the two dimensional distribution of the local stretching rate 
in a 2D surface flow with compressibility C = 0.21. The spectrum of the scaling exponents of 
the intermittent distribution of the coarse grained mass at small scales has been numerically 
computed. Our results indicate that the mass moments saturate at order ric — 0.86 in qualita- 
tive agreement with the prediction based on a Gaussian, white in time model of compressible 
flow. 

Direct measurement of the coarse grained mass density at small scales confirms the results 
within the error bars up to order n ~ 1.3. The events of high mass concentration at small 
scales have been identifled by calculating the probability density of the maximum number 
of particles on grid resolution. The multi-fractal spectrum obtained from the coarse grained 
moments for those specific events is compatible with the existence of a saturation exponent 
as predicted by the stretching rate analysis. 

* * * 
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